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O ! Various aspects of the Haldane-Shastry spin chain with 1/r exchange, and its various 



generahzations, are reviewed, with emphasis on its Yangian quantum group structure. 



O ' cind the interpretation of the model as the generalization of an ideal gas (of "spinons" ) to 

^ . the case of fractional statistics. Some recent results on its dynamical correlation function 
are discussed, and conjectured extensions of these remarkably simple results to the SU{n) 

^ I model, and to the related Calogero-Sutherland model with integer coupling constant, are 

Vh I presented. 
^' 

Q ! 1. Introduction 



Most of our understanding of quantum many-body physics is based on solvable models 



;_j ■ - the ideal Bose and Fermi gases, and harmonic oscillators. These are the paradigms, 

^ solvable in full detail, and treatments of more complicated, interacting systems usually 

aim to find a description of the system that is close to one of these paradigms. There 
are very few other fully- solvable models; the only general class of non-trivial solvable 
quantum models are the integrable one-dimensional models that can be traced back to 
the S = 1/2 Heisenberg chain solved by Bethe |I| in 1931. Starting with Bethe, there has 
been success in calculating the energy eigenvalues and thermodynamics of these models, 
and the algebraic structures that make them solvable have been identified as the quantum 
Yang-Baxter equation and "quantum groups" , but progress in explicit calculation of their 
correlation functions has so far been very limited. 

Bethe's model and its generalizations generally involve contact or delta-function inter- 
actions. A different class of integrable models based on inverse- square interactions was 
introduced around 1970 by Calogero p[ and Sutherland H. However, this class of models 



involved only gases of spinless impenetrable particles, and did not have the richness and 
lattice generalizations of the Bethe family. On the other hand, certain (static) correlations 
functions were explicitly found 0. A few years ago, spin-chain relatives of these models 
were discovered (the "Haldane-Shastry" model [§,0), and it was subsequently found how 
to put internal spin degrees of freedom [§|, ^ |], 0] into the Calogero-Sutherland model 
(GSM). 

In this paper, I will focus on the inverse-square-exchange spin-chain, and its trigono- 
metric (1/sin^) variants @, |^, mentioning in passing also its hyperbolic (1/sinh ) and 
elliptic (1/sn^) variants [|10|, which provide a continuous interpolation to Bethe's model. 
The trigonometric models are remarkable because they are explicitly solvable in much 
greater detail than Bethe's model, and share many of the characteristics of the ideal gas, 
but with non-standard or fractional statistics, which can occur in spatial dimensions be- 
low three. In particular, it seems that the action of the physically-relevant local operators 
(such as the spin operator on a given site of the spin chain) on their ground states excites 
only a finite number of elementary excitations. This parallels the action of a one-particle 
operator on an ideal gas ground state, and greatly simphfies the calculation of corre- 
lation functions in terms of the "form factors", the matrix elements of local operators 
between the ground state of the system and eigenstates characterized by finite numbers 
of elementary excitations. 

In the rational limit (pure 1/r^ interactions), it has proved possible to obtain simple 
explicit closed-form expressions for the thermodynamic potentials ||TT| , |12|, and very re- 
cently, for ground-state dynamical correlation functions |T^, which are again simple, but 
non-trivial. I believe that these models will finally be solved in full detail, approaching 
that with which the ideal gas can be solved, including, for example, full correlation func- 
tions at all finite temperatures. This should provide the first fully-developed extension to 
the standard paradigms. 

In this presentation. Section (2) will provide a self-contained introduction to the Yan- 
gian "quantum group" [jl4|, |15|, |16] symmetries of these models; I will mainly restrict the 
discussion to the S = 1/2 spin-chain models with SU{2) symmetry, and try to avoid the 
more esoteric mathematical characterizations (Hopf algebras, coproducts, etc.), disguising 
them in a more pedestrian physicists' terminology. Section (3) describes the application 
of this to the eigenspectrum of the trigonometric models, and Section (4) describes the 
recent remarkably simple results |TB| for some dynamical correlations of the 5* = 1/2 chain 



(and the A = 2 GSM), and presents new conjectured (and certainly correct) generaliza- 
tions to the SU{n) version of the chain, and to the GSM with arbitrary integer coupling 
constant A. 



2. The Yangian "Quantum Group" and Integrable Heisenberg Chains 

In this section, I will present a self-contained outline of the Yangian "quantum group" al- 
gebra, and its application to the Haldane-Shastry model (HSM) and its variants, pointing 
out a number of open questions. The HSM contains a novel realization of the Yangian 



algebra, originally characterized by Drinfeld |Tj] in the context of the algebraic Bethe 
Ansatz [|1^], but a mathematical structure in its own right. I will not give a compre- 
hensive history of this elegant mathematical edifice; a concise account by Kirillov and 



Reshetikin |T^ references many of the original works. An account of the representation 
theory of the Yangian is given by Chari and Pressley |19 



The original integrable model is the S = 1/2 Heisenberg chain, solved by Bethe in 
1931 i: 

H = J2P,,^,, P^^ = l + 2S,-S,. (2.1) 

i 

The integrability of Bethe's model derives from an underlying "quantum group" algebra 
called the Y{gl2) Yangian fl^: let T(m) be a 2 x 2 matrix with non-commuting operator- 
valued entries that depend on a "spectral parameter" u, act in the Hilbert space of the 
spin chain, and obey the algebra 

Ri^TiT^ = T'TlRl^ (2.2) 

where T^ is the 4 x 4 operator- valued matrix {T{ui)^l), T^ is (1®T(m2)), and R^^(-Ui,M2) 
is a 4 X 4 c-number matrix defined in the direct product space V^ ® V^ of two 2x2 
c-number matrices. Consistency of the algebra requires that R satisfies the quantum 
Yang-Baxter equation (QYBE) R^^p^^^p^^s ^ r23k13r^i2 -^^ Vi(g)V2 V3. The Yangian 
Y{gl2) corresponds to the rational QYBE solution 

R^"^ = {ui-U2)l + hP^^; (2.3) 

here P"*^^ is the exchange matrix where (for c-number matrices) P^^(A®B)P-'^^ = (B® A), 
and h is the quantum parameter, which in this case sets the scale of the spectral parameter. 
The Yangian algebra can be written as the commutation relation 

{u - v) [T''^{u),T^\v)] = h{T'<^{v)T^\u) - T^^(u)T"^(i;)) (2.4) 

It is commutative in the "classical" limit /i — > 0. The algebra implies that 

[t{u),t{v)] = Q, t{u)=i:i[T{u)l (2.5) 

so the t{u) are an infinite set of commuting operators. Another consequence is that 
[T(m), DetQ|T(f)|] = , where DetQ|T(t')| is the quantum determinant pO[| 

DetQ|T(M)| = T^\u)T^\u - h) - T^\u)T^\u - h). (2.6) 

The quantum determinant commutes with all elements of the algebra and is analogous to 
a Casimir operator. In the "classical limit" /i — >■ 0, the quantum determinant reduces to 
the usual determinant of c-number-matrices. 

It is consistent to impose the asymptotic condition T(m) -^ 1 for n ^ 00. An asymp- 
totic expansion can then be defined: 



where <t are Pauli matrices, and {a^} are the commuting generators of an infinite set of 
Abehan subalgebras of y ((7/2): [flm,o„] = [Jm,o.n] = [J^,an] = 0; they alone determine 
the quantum determinant: 

DetQ|T(u)| =0(m)0(m-/i). (2.8) 

The infinite-dimensional non-Abelian 1^(5/2) Yangian subalgebra, is completely generated 
by Jo, the generator of the s/2 Lie algebra (of SU{2) generators, so [Jq, J^] = ie"-'"^J^) and 
the additional generator Ji. All other operators can be expressed in terms of these: for 
example, 

J° = i/iJo-/o; J° = /iJo-Ji; ... , (2.9a) 

J2 = -iJixJi + lh\Jo-Jo)Jo; ... . (2.96) 

Note that, as a consequence of the Yang-Baxter relation, [J^, J°] = 0; also [Jq, J°] = 0, 
but [Ji, jO] ^ 0. 

The requirement that T(n) obeys the Y{gl2) algebra imposes consistency conditions 
on the Y{sl2) generators; the first non-trivial condition is the Serre relation 

[JlJl] = [r,,Jl]. (2.10) 

In fact, this relation is sufficient to ensure that Ji generates 1^(5/2), and can be used to 
compute the value of /i^; either sign oi h = ii/Zi^ can be consistently chosen. If the 
"classical limit" h ^ is taken, the J„ obey the infinite-dimensional Lie algebra 3/2+: 

[JLJn]=^^''''J:;.+n, m,n>0. (2.11) 

(A familiar example of this algebra is the subalgebra of any SU{2) Kac-Moody algebra 
defined by its non-negative modes). 

The fundamental representation oiY {gl2) (called the "evaluation" in the mathematical 
hterature [0, |ll) is 

Ti(m) = 1 + i; Pi = U + ^-5i (2.12) 

u 

where 5*1 is a S* = 1/2 spin. In this irreducible representation, Ji = 0, and the quantum 
determinant is {u + h)/u. 

If T(ti) is a representation of Y{gl2)., so is f{u)T{u — a), where f{u) is a c-number 
function {f{u) —>■ 1 as u —>■ 00) and a is a spectral parameter shift. The corresponding 
change in 1^(^/2) is Ji ^ Ji + aJo- 

Like Lie algebras, "quantum group" algebras have a fundamental property that allows 
larger representations to be constructed from tensor-products of smaller representations 
(the "coproduct"); however, unhke Lie algebras, the tensor-product operation is "non- 
cocommutative", which means that the result of a tensor product depends on the order 
in which it is carried out. In the case of Y{gl2), if Ti{u) is a representation acting on 
a Hilbert space Tii and T2{u) is another representation acting on another Hilbert space 
H2, then the matrix product 

T(n) = T^iu)T2iu) (2.13) 



is also a representation, acting on a Hilbert space Tii (8) 7^2, and 

DetqlTHI = (DetQ|Ti(M)|)(DetQ|r2(M)|). (2.14) 

For the ^(5/2) subalgebra, the tensor product representation is 

Jo = Jo(l) + Jo (2) (2.15a) 

(the usual "law of addition of angular momentum" ) and 

Ji = Ji(l) + Ji(2) + ihJoil) X Jo(2). (2.156) 

Note that reversing the ordering of Jo(l) and Jo (2) defines an alternative tensor prod- 
uct with h ^ —h; either definition can be used, which is why only /i^ is fixed by the 
Serre relation. In Y{gl2) the two alternative tensor products are right and left matrix 
multiplication. 

A (Yangian) highest-weight state (YHWS), is a state that, for all u, is annihilated 
by T^^(u), and is an eigenstate of T^^{u) and T'^'^iu) with eigenvalues 4>i{u) and 4>2{u), 
where 

(DetQ|T(u)|)|YHWS) = 0^(n)0+(n - /i)|YHWS). (2.16) 

Every finite-dimensional representation has at least one highest-weight state. It is irre- 
ducible ii there is only one YHWS; conversely, an irreducible representation can be derived 
from each YHWS contained in a reducible representation. To each such YHWS can be 
associated an invariant subspace of all the states that can be generated from it by suc- 
cessive action of elements of the algebra. In contrast to what happens in the unitary 
representation theory of Lie groups and their associated Lie algebras, the representation 
of Y{gl2) within this invariant subspace may still be reducible, as the subspace may con- 
tain other YHWS (and their invariant subspaces) which are eigenstates of the quantum 
determinant with the same eigenvalue. 

An irreducible representation containing (only) a given YHWS is obtained by con- 
structing its invariant subspace, and then projecting out the invariant subspaces of all 
other YHWS contained within it. This is not necessary in the important special case 
of full reducibility, associated with "quantum symmetries" , when the invariant subspaces 
associated with each independent YHWS are always irreducible, and orthogonal to each 
other. Reducibility without /u// reducibility, where the YHWS |2) is contained within the 
invariant subspace generated by YHWS |1), means that there are non- vanishing matrix 
elements of the form (2|/(T(u))|l), but all matrix elements of the form (1|/(T(m))|2) 
vanish. 

The representation theory of Y{sl2) has been described in detail in [jl9|. The basic 



theorem is that all finite-dimensional irreducible representations are isomorphic to an 
irreducible representation derived from the maximal- Jq YHWS of some tensor-product 
of fundamental representations. This means that the eigenvalues of T^^{u) and T'^^{u) 
acting on a YHWS which generates a finite-dimensional invariant subspace have the form 

0+(w) = /(n)P+(w + h); 0+(n) = /(M)P+(n), (2.17) 



where P^{u) is a finite-dimensional polynomial called the Drinfeld polynomial. Let 

Nn /n-1 



,(n) 
\a — C 

n>li=l \u=0 



PHu)=l[Il[U[^-^i+^h]\, (2.18 



be the unique decomposition of the roots of P^{u) into "strings" of consecutive roots 
with spacing /i, obtained by successively factoring out strings, starting with the longest 
ones present (this fixes how multiple roots are treated). The dimension and SU{2)- 
representation content of the irreducible representation of Y{gl2) derived from the YHWS 
can immediately be identified from this form: it is equivalent to a tensor product of SU{2) 
representations where S = n/2 occurs Nn times, with a total dimension (2^^) (3^^) . . . ((A;+ 
1)^'=), where the YHWS has Jg = k/2 = J2n'^^n/2. If only 1-strings are present in 
the Drinfeld polynomial of a YHWS, the invariant subspace generated from it must be 
irreducible. 

The Algebraic Bethe Ansatz |T^ is the solution to the problem of finding the eigen- 
vectors of the commuting operators t{u) within an irreducible finite-dimensional repre- 
sentation of Y{gl2): the eigenstates with eigenvalue ty{u) which are SU{2) highest-weight 
states with Jg = {N/2) — M^, are given by 

|z.)=(^nTi2(t;f))j|YHWS) (2.19) 

where the M^'^^ "rapidities" {fi*^ } are solutions of the Bethe Ansatz equations (BAE) 

4>t{u)Qu{u ~h) + 0+(m)Q^(m + h)= t^{u)Q^{u) (2.20) 

and Qu{u) is the polynomial 

QAn) = Ei^-^h- (2-21) 

The eigenstates of t{u) form an orthogonal basis only if [(Jg • Ji), (Jo " <^i)^] = 0. 

The Bethe model on a periodic chain of N sites with 5=1/2 spins corresponds to a 
Y{gl2) representation that is a simple tensor product of fundamental representations with 
no spectral parameter shifts: 

T(m) = Ti(m)T2(m)...T;v(w). (2.22) 

It is irreducible, with (piiu) = {{u + h)/u)^ and 02^ (m) = 1. The limit as u ^ of u'^t{u) 
is h'^ exp iK, where exp iK is the lattice translation operator, and we may write 



Tiu) = T^iu)T2iu)...T^{u). 

is irreducible, with 4>t{u) = {{u + h)/u)^ and 02^ (m) = 1. The limit as u - 
h^ exp iK., where exp iK is the lattice translation operator, and we may 



«^ 



t{u) = h^e'^ exp Y. w"^n+i , (2-23) 



vn.>l 



where {Hn} are a set of local Hamiltonians. Then [Hm, Hn] = [Hm, Jo] = [Hm, exp iK] = 
0, and H oc H2. The 1^(^/2) generator may be written 

J^ = ihY,SiX Sj. (2.24) 

i<j 



In the conventional Bethe parameterization, u is the rapidity parameter if the choice h = 
i is made; in this case Ji is Hermitian. 

For finite N, [expiK, Ji] ^ 0, and [i^m; ^i] 7^ 0. However, in the thermodynamic 
limit N —^ oo, the model acquires a full non-Abelian quantum, symmetry: the Yangian 
generator Ji now commutes with the lattice translation operator, and with all the Hn- The 
infinite-dimensional representation of the Yangian becomes fully reducible into orthogonal 
subspaces in which the Hn are diagonal. Furthermore, the antiferromagnetic ground 
state becomes the unique state that is singlet under the action of the Yangian. In the 
thermodynamic limit a spectral parameter "boost" operator B [^, where [Ji,i?] = hJo, 
is given by 

^ = E(^ + l)Pii+i- (2-25) 

i 

This operator, which is incompatible with periodic boundary conditions, also generates the 
family of commuting Hamiltonians with quantum symmetry from the lattice translation 
operator: [B,expiK] = hH2expiK, [B,Hn] = nhHn+i- 

The most general representation of the Y{sl2) Yangian based on a simple tensor prod- 
uct of fundamental representations (including spectral parameter shifts) is 

^1 = E 7.5^. + I E ^*. S^ X S, (2.26) 

where Wij = ihsgn{i — j), and the {7,} are arbitrary. It is an instructive exercise to 
substitute a form of this type into the Serre relation, without making any assumptions 



about Wij, except that it can generally be taken to be odd. The result [|T5| is that this is 
a representation of the Yangian generator if 

WijWjk + WjkWki + WkiWij = h"^, i ^ j ^ k. (2.27) 

The general solution of this is 

^.^. = ^h (^1±^] (2.28) 

V ^i ^j I 

where {2:4} are a set of A^ arbitrary, but distinct, complex parameters. The simple tensor- 
product representation is recovered in a limit where \zil Zj\ — > 00 for i > j. 

Note that 1^(^/2) representations with an invariant subspace where groups of 25* con- 
secutive 5* = 1/2 spins remain in the symmetric spin-S" state can be obtained from the 
simple tensor product by suitably choosing the associated spectral parameter shifts 7j to 
be a 2S'-string, and this is the basis for the extension of Bethe's model to chains of spins 
with S > 1/2 [Q. However, this requires the stronger property {wij^ + h"^ = 0, which 
eliminates more general representations based on (2.26) and (2.28) for such spins. 

The 1^(5^2) representation based on (2.28) is evidently derived from a more general 
class of representations of Y{gl2) than those obtained by a simple tensor product. Instead, 
following Ref.[|T^], they can be obtained as follows. Let {74} and {zi} be variables that 
commute with each other, and with the spin variables, and let Kij be a permutation 



operator that permutes their labels: Kij'ji = 'jjKij, etc. These permutation operators do 
not affect the spins: [K^, S^] =0. Now consider the quantities 



These variables commute, and obey a "degenerate affine Hecke algebra" ||16|| : 

Kii+i^i - 7i+i^ii+i = h; [% 7j] = 0. (2.30) 

(The role of Hecke algebras in this context has been stressed in Ref. [0.) This implies that 

[K,j,A{u)] = 0, A{u) = l[{u-%), (2.31) 

i 

and hence that A(m) is a c-number function of u and the {7^, Zi}. Since the 7^ commute 
with each other and with the spins, a representation of Y{gl2) is clearly given by 

T{u) = Ti(m - 7i) . . . T^iu - fjv), (2.32) 

where A{u)T{u) is a polynomial of degree A^ in u. 

Now let n"*" be a projection operator into the subspace that is fully symmetric under 
simultaneous permutation of spins and the parameters {7j,2;j}, so Kijll~^ = PijU^. The 
Hecke algebra guarantees the property 

H+T(m)H+ = T(m)H+. (2.33) 

This ensures that the totally-symmetric subspace is an invariant subspace of T(m), and 
hence that the projection of (2.32) into this subspace is also a representation of Y{sl2)- 
The representation can now be evaluated |]16| within the totally-symmetric subspace as 



T{u) = 1 + /i E Pi((^ - L)~%, (2.34a) 

where Lij are elements of a A^ x A^ "quantum Lax matrix" 

Lij = lAj + h{l - 5ij){l - {zi/z,))-^Pi,. (2.346) 

The 2^ states of this representation are all eigenstates of the quantum determinant with 
the same eigenvalue A{u + h)/A{u) where A{u) is the polynomial of degree N given by 
det |m — L°| and L^^ is the eigenvalue of Lij acting on the fully-spin-symmetric YHWS 
with J§ = N/2 and Pij = 1. This YHWS has the Drinfeld polynomial P+{u) = A{u). 

For generic values of the {74} and {zi}, the representation is irreducible. We can 
now ask whether for some special choices of these parameters, the representation becomes 
fully-reducible as a direct sum of smaller representations, which would imply that there is 
a non-trivial Hermitian operator H that commutes with both Jq and Ji. We can postulate 

8 



such an operator to be of the form H= J2i<j hijPij, and look for solutions of the condition 
[H, Ji] = 0. This is satisfied provided 

7i - Ij = 0; hij oc [wfj + h^)] ^ Wijhij = 0. (2.35) 

We also require that H is Hermitian (hij real). We find [^ two families of translationally- 
invariant Hamiltonians with quantum symmetry, where [expiK, Ji] = [expiK, H] = 0; the 
Hamiltonian is even parity under spatial reflection, while Ji is odd-parity, with ( Ji)^ ex Ji, 
which ensures that any reducibility is full reducibility. The solutions are 7j = 0, and hij 
oc l/d{i — iY, where in the hyperbolic models, defined on an infinite chain. 



w. 



«j 



i/i coth(/t(i — j)); d(i) = K ^ sinh(fi:j) (2.36a) 



with real Zi+i/zi = exp2K, which can be either positive (k real) or negative [k — iTT/2 
real). In the limit k -^ oo, the hyperbolic model becomes the N = oo Bethe model; in 
the limit k — *> 0, and ih = 2k, so h ^ 0, we get a "classical limit " where Y{sl2) becomes 
the infinite-dimensional Lie algebra 5/2+, and we obtain the rational model 

w., = {t-j)-'; d{j)=j. (2.366) 

The other family of solutions are the trigonometric models, obtained from the hyper- 
bolic models by the replacement k, -^ in/N, with 

Wij = hcot{7T{i - j)/N); d{j) = (N/n) sm{7f j/N), (2.37c) 

defined on a finite periodic chain on A^ sites, with d{j) being the chord distance between 
lattice sites equally-spaced on a circle. For the choice of real h, the A^ roots of A{u) form 
an A^-string along the real axis, centered on the origin. The trigonometric model then has 
{LijY = Lji, and hence {T'"^{u)y = T^"(u*). 

In all these models, H oc H2 is the first member of a family of commuting constants 
of the motion {H2,H3,H4, . . .} that commute with Jq, Ji, and expiK. The eigenstates 
of these operators form irreducible representations of Y{sl2)- For example |10[| , 

H^oc Yl —^ Si ■ Sj ^ Sk, Zij = Zi- Zj, (2.37) 

but the systematic construction beyond i/4 |]15| has not yet been elucidated. So far, no 



generalization of the spectral-parameter "boost" operator from the Bethe limit to the 
hyperbolic and trigonometric models has been found. The origin of the Hamiltonian 
constants of the motion in the models with quantum symmetry is conceptually rather 
different from that in the finite- A^ Bethe model. They must commute with all elements 
of T{u), but the only elements of Y{gl2) with this property derive from the quantum 
determinant, which is a trivial c-number in these representations. Thus the Hamiltonians 
with quantum symmetry cannot be expressed as functions ofT{u), and are independent 
objects. 



The Yangian symmetry of the hyperbolic models is essentially similar to that of the 
Bethe model: it occurs in the thermodynamic limit, and Ji is made Hermitian by the 
usual BAE choice h = i. On the other hand, the realization of Yangian symmetry in the 
trigonometric models is essentially new, as it occurs in a finite periodic system, which 
has discrete, highly-degenerate energy levels that can be classified as finite-dimensional 
irreducible representations of 1^(5/2) . In this case Ji is made Hermitian by choosing h 
real. The two variants are only connected via the "classical" rational model with h = 0. 

The consequences of Yangian quantum symmetry are very different in the two cases; 
there is an analogy to the difference between the non-compact Lorentz group in (1+1)- 
dimensions and SU{2), which have symmetry algebras corresponding to different real 
forms of the same complex s/2 Lie algebra. In the hyperbolic case, the representation 
theory involves the Bethe "rapidity strings" with spacing h = i in the imaginary rapidity 
direction, which are related to bound states. In the trigonometric case, the string spacing 
is in the real direction, and is related to the quantization of the momentum of free particles 
with periodic boundary conditions in units of 27c/N, and to a generalization of the Fock 
space structure of such a system. 

The trigonometric models, and their rational limit, were independently discovered 
by Haldane Q] and Shastry [0. The extension to the hyperbolic model, was found by 
Inozemtsev [0, who also proposed the integrability of an elliptic variant with hij oc 



p{i — j), where p{u) is the doubly-periodic Weierstrass function with periods A^ and 
in/K. This is a version of the hyperbolic model, made periodic on a finite chain of A^ 
sites, and (like the finite-A^ Bethe model (k -^ 00)) does not have quantum symmetry. 
A proof of integrability is missing, but two odd-parity operators that commute with the 
Hamiltonian are given in |T^, making integrability very plausible. One linear combination 



of these corresponds to H^ in the trigonometric and hyperbolic limits, the other, which 
corresponds to Jo ■ Ji, is (where ({u) is the elliptic zeta function) 

E (C(^ - j) + C(J -k) + Cik - i)) S, X S, ■ Sk. (2.38) 

I have not been able to find a Y{sl2) generator Ji where Ji ■ Jq gives (2.38). I therefore 
conjecture that the integrability of this model may involve a yet more general "quantum 
group" algebra with two "quantum parameters" h = 2m and h' = 27: /N (which would 
be scale parameters for two distinct spectral parameters), which only degenerates to the 
Yangian when either of them vanishes or becomes infinite (the Bethe limit). This would 
correspond to a "double quantization" of 5/2+, or a further "quantization" of Y{sl2)- 



Indeed, Inozemtsev [^ suggests that two spectral parameters play a role in this model. 
A better understanding of the elliptic model is clearly needed. 

To conclude the formal discussion of algebraic aspects of these models, I note that. 



as in the case of Bethe's model [^, there is a straightforward extension from the Y{gl2) 
quantum group with SU{2) symmetry, to the Y{gln) quantum group, where T(m) is an nx 
n matrix, with SU{n) symmetry. A further extension is to the graded or supersymmetric 
generalization, Y{glm\n), with SU{m\n) supersymmetry, where a site can be in one of m 
states with even fermion number or n states with odd fermion number. Let cL create on 
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site i a particle of species a, which may be fermionic or bosonic, and impose the constraint 

^cLqc, = 1, alH. (2.39) 

a 

The exchange operator becomes 

Pij = H 4aC]pCif3Cja, (2.40) 

a/3 

and the generators of the slm\n graded Lie algebra and its Yangian extension can be given 
as the traceless operator-valued matrices 

Jf = E (4q, - (m + n)-'6^^) , (2.41a) 



Ji = I E Wijcl^c]^Ci^Cjf3. (2.416) 

In particular, the graded trigonometric exchange model with {m\n) = (1|2) corresponds 
to the "supersymmetric t-J model" variant found by Kuramoto 



3. Spectrum of the Trigonometric Haldane-Shastry Model, and its Interpre- 
tation as a Generalized Ideal Gas 

In this section I will focus on the remarkable properties of the trigonometric model, 
with its realization of Yangian quantum symmetry in a form compatible with periodic 
boundary conditions. The Hamiltonian, normalized so the spin-wave velocity is Vg in 
units where the lattice spacing a = h = 1, is 

N^ ^^sm\nit -j)/Ny ^ ' ' 

Let us first consider the case where the symmetry group is SU {1\1), when the model is 
equivalent to a lattice gas of free spinless fermions, with creation operators a] = Ci2Cii, so 

Pij = 1 — (al — aj){ai — aj) = — 1 + (oj — aj){al — Oj). (3.2) 

The energy levels of this model are characterized by a sequence of A^ — 1 Bloch-orbital 
occupation numbers {rii, n2, . . . , um^i}, taking values or 1, so 

N-l N-1 

E=J2 Unm - I); e'"" = n e*"'""™ (3-3) 

m=l m=l 

where 



,'27TV, 
Cm 



[^-j^jmim-N); k^ = 27im/N. (3.4) 
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All eigenstates have a two-fold degeneracy, corresponding to the two possible occupation 
states of the translationally-invariant Bloch orbital, which has zero energy. This is the 
supersymmetry: all states are s/i|i doublets, and there is no non-trivial Yangian extension 
in this case. The ground state of H is the state with maximum occupancy, {11^} = 
{111... 11}, and the ground state of —H has minimum occupancy {000... 00}; The 
replacement 1 <->■ in any configuration maps E to —E and K to 7r(A^ — 1) — i^. 

It is a remarkable fact that this set of energy-momentum levels is the complete set of 
levels contained in the general SU{'m\n) exchange model, but in the general case they 
have different and much larger Yangian multiplicities. Specializing to the SU{n\Q) and 
S'f/(0|ra) cases, some of these multiplicities are zero, and the energy level is absent. The 
selection rule for allowed multiplets in the case SU{n\{)) is that occupation patterns con- 
taining a sequence of n or more consecutive I's are forbidden; in the SU{Q\n) case, n or 
more consecutive O's are forbidden. Each allowed sequence corresponds to an irreducible 
Yangian multiplet. In general, the spectrum of H with SU{m\n) symmetry is the same 
as that of —H with SU{n\m) symmetry. 

The above rule was found empirically by examination of numerical diagonalization 
results [|15|, but I now specialize to the simplest case, SU{2\Q) = SU{2), where the rule is 
that occupation patterns with consecutive I's are forbidden, and make contact with the 
Y{sl2) representation theory |[19[ . 



Recall that each YHWS labeled by z/ has an associated Drinfeld polynomial, P^{u), 
and that A{u)T{u) is a polynomial, with quantum determinant A{u — h)A{u + h). 
Thus there is a polynomial fu{u) where 4>i{u), (the eigenvalue of T^^(m)), has the form 
Uu)P+iu + h), 0+(m) = Uu)P+{u), and 

fu{u)Uu - h)P^{u + h)P^{u -h) = A(n + h)A{u - h). (3.5) 

Since the roots of A(u) are an A^-string, it is easily seen that P^{u) must be a factor of 
A(m), so A(m) = P^{u)g,y{u), where gu{u) is a polynomial with roots at the roots of A(u) 
not contained in P+(u), and 

fy{u)fu{u + h) = g„(u - h)g^(u + h). (3.6a) 

This is a polynomial equation with the elementary solution 

g{u) = {u — a){u — a + h), f{u) = [u — a — h){u — a + h) . (3.66) 

Any product of such solutions is a solution. This shows that gviu) must be a product of 
2-strings. 

We now have the recipe |T^ for constructing the possible Drinfeld polynomials (in 



fact, there is one YHWS in the spectrum of the trigonometric model corresponding to 
each allowed Drinfeld polynomial [l^). First partition the A^-string of roots of A(u) into 
1-strings and 2-strings. Between each consecutive pair of roots place a 1 if they belong 
to the same 2-string, and otherwise. This gives a binary sequence of length A^ — 1, 
with the constraint that there are no consecutive 1 's, as found empirically. The 1-strings 
are the roots of the Drinfeld polynomial; if an extra is added at each end of the binary 
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sequence, the 1-strings are located between each consecutive pair of O's. From the Yangian 



representation theory, we conclude (in agreement with the empirical findings ||T5[) that a 
sequence of M + 1 consecutive O's represents an independent S = M/2 degree of freedom. 
There is a simple physical interpretation: each root of the Drinfeld polynomial repre- 
sents the presence of a 5* = 1/2 "spinon" excitation. A m-string of such roots represents 
m spinous "in the same orbital" with a rule that their spin state must be totally symmet- 
ric. If the total number of spinous present is Ngp, there are M = {N — Nsp)/2 I's in the 
occupation pattern, which serve to separate the A'^orb = M + 1 "orbitals" into which the 
spinous are distributed. The change in the number of available "orbitals" as the spinon 
number is changed obeys ANorh/^Ngp =—1/2, which allows spinous to be identified as 



excitations with "semionic" fractional statistics, in between Bose and Fermi statistics ||25 

If {Jj} are the positions of the M non-zero entries in the length-(A^— 1) binary sequence 
{rij}, "spinon orbital occupations" {rii^}, 1 < i < M + 1, a = ±1/2, are given by Jj+i — Jj 
= 2 + riji + riji, with Iq = —1, and /j\/+i = A^ + 1. We can then treat the {riia-} as 
independent variables subject only to the constraint 

N,p = J2n^a = N-2M. (3.7) 

Each distinct configuration {rii^} satisfying this constraint corresponds to an eigenstate 
of the system. The energy and momentum can now be expressed in terms of the "spinon 
orbital occupations": it is convenient to relabel Ui^, with 1 < i < M -|- 1, as n^a, where k 
is a crystal momentum in the range —ko<k<ko,ko = ttM/N, then if 

e{k) = ^{ki-ey, v{k) = v,{ko~\k\y, (3.8) 

vr 
the spectrum is given by 

E = Emn + Yl ^ik)nka + ^ Yj ^(^ ~ k')nk^nk'a', (3.9a) 

ka ka-,k'<j' 



e 



iK 



_^)M-Qg.fcn,.. jz^J2^^^^_ ^3 9^) 



ka k<j 



Here Emn = TTVsi2N{N^ - 1) + 4:M{M^ + 2) - 3MN^)/QN^. From this spectrum, it 
is straightforward to derive the thermodynamic potentials |]ll| in the limit N ^ oo, 
including a Zeeman coupling —hJ^. 

The spinon orbitals can be parameterized by a rapidity x in the range — 1 < x < 1, 
(where the spinous in that orbital have velocity Vgx), and have mean occupation numbers 

Uaix) = exp —P[e{x) — crh{l + fi{x))], (3.10) 

where e{x) = {7iVs/4:){l — x^) and 

smh{phfx{x)/2) 



smh{ph/2) 
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exp(-/3e(x)). (3.11) 



This gives an easily-solved quadratic equation for exp{j3hfi{x)/2). The free energy per 
site is given by 



-m..)4£^-'- f°"S;^'"^ - 



(3.12) 



In the absence of a magnetic field, the entropy per site is 

s{(3,h = 0) = kB f dx{ln[2cosh/3e{x)]-f3e{x)tanhf36{x)), (3.13) 

which unexpectedly is even in (3, so is the same for the ferromagnetic and antiferromagnetic 
models. 

The identification of spinous as semions is supported when the wavefunctions of the 
trigonometric SU{2) model are examined. A class of polynomial wavefunctions of the 
Calogero- Sutherland model of a non-relativistic gas with 1/ sin^ interactions was discov- 



ered [§|, g to also give a class of "fully-spin-polarized spinon gas" [^ wavefunctions of 
the trigonometric model, and these are now identified with the YHWS states. If Zi (with 
{Zi)^ = 1) are the complex coordinates of the M sites with ctj = —1/2, the wavefunctions 
have the form 

^.{{Z,}) = <|.,({Z,,})vI/o({Z,}), (3.14) 

where $j,({Zj}) is a symmetric polynomial with degree Ngp in each Zj, and '^q{{Zj}) is 

M{Z^^) = I[iZ,-Z,fl[Z,. (3.15) 

The is essentially the same as the m=2 Laughlin fractional quantum Hall effect state for 
bosons p6|, and for N^p = 0, this is the unique Yangian singlet state. The symmetric 



poljTiomials $y are solutions of the eigenvalue equation 

(e-I^ + AE (^) 4) *'^'({..}) ^ .*'^'({..}). (3.16) 

with A = 2. (The solutions are known in the Mathematical Literature as the Jack Poly- 
nomials |2^). If {/j} are the positions of the non-zero entries in the binary "occupation 
number" sequence, the Taylor series expansion of the corresponding YHWS wavefunction 
has the form 

^=Y: C{{m,}) Y: (EiZpi^)^] , {^^} < U^} (3.17) 

{mj P \ i / 

where P{i) is a permutation, and where rrij < rrii^i, and {rrii} < {m'^} means that {rrii} 
can be reached from {m'^} through a sequence of pairwise "squeezing" operations rrii ~^ 
rrii + 1, rrij — ^ rrij — 1, with rrii < itt-j — 1, and rrik unchanged for k ^ i,j. 

The Nsp = Yangian singlet state may also be written in terms of the azimuthal 
spin variables ai: it occurs only for even N, and is the n = 2 case of the SU{n) singlet 
wavefunction where the Uj can take one of n ordered values {a}: 

¥^\{zi,a,}) =l[{z,- zjY^'^--^\rr-^--'^^^l[6{N,nN{a)), (3.18) 

i<j a 
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where N{a) = I]i 5(q;, ctj). Here the {zi} are a complete set of the A^'th roots of unity. 
Note that this state is essentially the same as the wavefunction for a filled Landau level 
of SU{n) fermions, and is hence explicitly 5'f/(ra)-singlet for arbitrary {zi}. 

The Ngp = 1 states are also particularly simple: they occur only for odd A^ = 2M + 1, 
and are generated by 

<if{z; {Z,}) = Yliz - Zi) \{{Z, - Z^f n Z^. (3.19) 

i i<j i 

Expanding this in powers of z gives a band of Bloch states with crystal momentum —tt/2 
< K + Mn < 7r/2, so the spinon band covers half a Brillouin zone. The localized spinon 
wavefunction "^Iz; {Zi}) is essentially the quasihole excitation of the m = 2 bosonic 
Laughlin state, and from this perspective, clearly describes a semionic excitation. If the 
parameter z is chosen to be a lattice coordinate Zi, the localized spinon wavefunction with 
ai = a can be rewritten in terms of spin variables as 

6{a,a,)¥;f\{z„af,jy^t}). (3.20) 

This shows that the spinon is completely localized on the lattice site, and induces no spin 
polarization of its local environment. While there are N such fully- localized spinon states, 
they are fundamentally non-orthogonal, since the expansion in orthogonal Bloch states 
shows there are only (A^ + l)/2 independent Nsp=l states. 

To end this section, I discuss the extension of the state-counting from SU{2) to the 
general SU{n\m) case. The Yangian counting rules for the SU{2) model are very simple. 



and allowed the thermodynamic functions to be explicitly obtained |Tl[] in closed-form in 
the thermodynamic limit, but become much more complicated in the general SU{m\n) 
case. Recently Sutherland and Shastry |12[ showed how to recover the SU{2) results for 



the thermodynamics, and generalize them to SU{m\n), from a strong-coupling limit of 
the exchange-generalization of the Calogero-Sutherland model [^, 0) §) 0- 

The spin chain degrees of freedom are essentially those of the spin-1/2 Fermi gas with 
the charge degrees of freedom removed. From this viewpoint, the spinon is the spin-1/2 
fermion with the charge degrees of freedom factored out, leaving "half a fermion", and 



hence a semion. I will interpret Sutherland and Shastry's result [|T2[ as showing how to 
"put back" charge degrees of freedom into the spin chain to recover a spectrum with the 
familiar degeneracies of the ideal gas with internal degrees of freedom. This facilitates 
the computation [|r^ of the thermodynamics in the general SU{m\n) case. 

Let bl, k = 1, . . . N — 1, he a set of harmonic oscillator creation operators, and add 
these degrees of freedom to the spin chain as follows: 



^iK' ^ ^iKs^iK, 



2_ /N-l \ 

Kc = ^[j:kbib,j, (3.216) 
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where expiKg is the spin-chain translation operator, and Pij is the SUlmln) exchange 
operator. 

Now let Uka = c'ka^ka be occupation numbers of a set of orbitals for particles of a 
bosonic or fermionic species a, with fc = 0, ±1, . . . ± oo, subject to the constraint that 

iV = ^n,„ (3.22) 

ka 

is fixed. From the result of Sutherland and Shastry, the spectrum of H' (with the con- 
straint (3.22)) is identical to that of 

H" = (S) EE 1^ - ^>fc«^fc'/3, (3-23a) 

^^^ ^ kk' af3 

K" = j^Y.knk., jf = Y.icL<^kp-{m + n)~'5-^). (3.236) 

ka k 

The spectrum of H" is periodic in K" -^ K" + 27r, corresponding to a shift rika — > ^fc+ia 
of the occupation number pattern. The precise statement is that the spectrum of H" in 
one period of K" coincides with that of H' . In the thermodynamic limit, the free energy 
of H' can thus be easily calculated, as can the free energy contribution from the extra 



oscillator modes; the difference is the spin-chain free energy |12]. However, this does not 
give a simple method for identifying the Yangian degeneracies of the discrete (finite A^) 
spectrum of H without the oscillator modes. The explicit expressions for the free energy 
are in fact only obtained in the limit N —>■ oo, when the Yangian algebra has degenerated 
to its "classical" s/min , limit. 



4. Dynamical Correlation functions 

Let us consider the state 5'./'|0), where |0) is the Yangian singlet ground state of the 
trigonometric chain. It is a spin-1 state, and is easily seen to be a linear combination of 
Yangian highest weight states with Ngp = 2. The action of S^ is to remove a down-spin 
coordinate at site i. Thus M = N/2 — 1 and 

M 

^{z,, {Z,}) = n(^^ - Z,f^oi{Z,}). (4.1) 

i=i 

The polynomial pref actor has degree 2, confirming that this state is composed purely of 
two-spinon eigenstates. Thus if we wish to compute the dynamical correlation function 

C{z-j,t-t') = {0\S-{t)S+{t')\0), (4.2) 

the only intermediate states that contribute to its spectral function are the Nsp=2 YHWS. 
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Recently, it has become possible to compute certain dynamical correlation functions 
of the Calogero-Sutherland model 0, ^ 

^ f P? , ^V A(A-l) 

with d{x) = (L/tc) sin{iTx/L), in the thermodynamic limit at fixed density p = N/L. The 
eigenf unctions of the CSM have the form 

^^'H{Z^})^j{{Z.}y, ^J = mZ^-Z,)'YliZ^y (A>0) (4.4) 

i<j i 

where {"^j} is the family of states that includes the ground state and Galilean boosts of it, 
Zi = exp 27iiXi/L, and ^'^^^({Zj}) is a symmetric polynomial solution of (3.16). While the 
apparent statistics can be modified with a singular gauge transformation, the "natural" 
statistics of this model are, from (4.4), evidently fractional, with statistical parameter 6 
= 7tX {i.e., particles carry charge 1 and flux ttA). The elementary excitations are particles 
with velocities greater than the speed of sound Vs = irhXp/m, and holes with velocities 
less than Vg- The holes carry charge — 1/A and flux —it, and have statistical parameter 6h 
= n/\. 

The calculations can be carried out at one of three special couplings (only two of which 
are non-trivial), and involve a mapping to the N ^ oo limit of a Gaussian dynamical 
N X N matrix-model, with either orthogonal (A = 1/2), unitary (A = 1), or symplectic 
(A = 2) symmetry. This reduces the problem to a complicated, but tractable Gaussian 
problem [P^. 

At integer couplings X = q, the natural particle statistics are Bose (even q) or Fermi 
(odd q). The ground state is ^j with J = —q{N — l)/2. The wavefunction of the state 
^(x)|J) is 

n(^ - Z,y<iIj{{Z,}), z = exp{2mx/L), (4.5) 

i 

which is composed only of eigenstates with just q hole excitations, so the spectral func- 
tion of the retarded single-particle Greens function will only contains contributions from 
intermediate states of that type. This correlation function will thus have the form 

{0\^^{x,t)^{0,0)\0) = p (n/%) l/.d^.DlV^^^-''^ (4.6a) 

P = ^mhVi, E = ^\mhv'l, mh = -m/q, (4.66) 

i i 

where fq{{vi\) is a form factor that must be calculated. 

Using the mapping to the symplectic matrix model, it has been possible [|13| to calcu- 
late this form factor for A = g = 2. The result is remarkably simple: 

l/.fa.^.)|--^( , ,'"'-f J . (4.7) 

8vs \{vi - vf){vi - vi) J 

17 



The significance of this in the context of the SU{2) trigonometric Haldane-Shastry chain is 
that the YHWS wavefunctions with Zi = exp{27rixi/L) are also eigenstates of the A = 2 
CSM, and the matrix elements of S^{t) between two YHWS are equivalent to those 
of '^{xi,t) between the corresponding CSM states. The CSM result can immediately 
be translated to the results for the spin chain in the rational limit, where the Yangian 
symmetry algebra becomes classical. 

While the calculation [|13| is complicated and indirect, the simplicity of the result 



suggests there should be a simple derivation, not based on the "accident" of a mapping 
to a Gaussian matrix model. It is very tempting to conjecture the extension of the result 
to the SU{n) chain, for which no such mapping is known. Examination of the action of 
the operator cl^Cip on the singlet ground state of the SU{n) chain leads to the conclusion 
that it produces a two-parameter family of YHWS states with one spinon (of "color" (3) 
moving with rapidity Xi, and a complex oi n — 1 spinous (with net "color" a) moving 
together (but not as a bound state) with rapidity X2- In the singlet ground state of the 
rational SU{n\0) chain state, the correlation function (0|X (t)XQ (0)|0), where X" = 

CiaCi/3, will thus have the form 



[-ly 



S'^s^yPCij, t) + i-i^r^,5^^(l - C(j, t)), (4.^ 



n n^ 



where 



Cim,t) = jf' dXl C dX2 |F„(X1,X2)|' (^e*(9(-l)— ^(-iW) (^^mx2)m-e(x2)t)Y ' ^ (4 g^) 

q{x) = —- e(x) = ^(l-x2). (4.96) 

n In 

Here F„(xi,X2) is the form factor that must be found. 

The asymptotic form of the correlations are straightforward to compute from bosoniza- 
tion (they are free-fermion correlations with the factor coming from charge degrees of free- 



dom divided out), or from conformal field theory |2^, and can be fit to a simple Ansatz 
based on (4.7). I therefore present the conjecture for the form factor of the rational 
SU{n\{)) chain (which is a rigorous result for n=2) : 

\{i - xl){i - xl) ) r(^)r(^) 



It has been verified |^0[ that this form fits the numerically computed static structure 
factor of the trigonometric SU{?)) chain with A^ < 18 extremely well, with very small 
finite-size corrections, leaving no doubt that this conjecture for ra > 2 is indeed correct. 

As another test of the reliability of generalizing (4.7) "by conjecture", based on its 
remarkable simplicity, the same type of arguments can be used to obtain the form factor 
for the retarded Greens function (4.6) of the CSM at a general integer coupling X = q. 
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In this case, the states contributing to the spectral function have q hole excitations with 
independent velocities. The resulting conjecture is 

\U{b^})? = l^B,[Jl{vl-vj)^ \\{{v.-v,f\ , (4.11a) 

where the multi-dimensional integral fixing the normalization Bg is hard to do. I have 
just learned of recent work by Forrester [^, who calculates the equal-time limit of this 
correlation function using quite different methods based on the celebrated Selberg trace 
formula. His expression is precisely the equal-time limit of the formula conjectured here! 
Moreover, Forrester's result provides the normalization as 

, ( r (i±A \ 




It is striking that the form factors are essentially Laughlin-type wavefunctions for anyons 
with the same statistics as the holes, now with the rapidities as coordinates. 

The remarkable property of the rational and trigonometric models is that the local 
operators such as the spin on a given site act on the ground state to produce only a very 
restricted class of excitations. There is a general selection rule, verified empirically ||32 



that the local spin operator Si acting on any eigenstate cannot change the spinon number 
Ngp by more than ±2. Such properties are very reminiscent of an ideal gas, and the most 
natural interpretation of the trigonometric and rational models is as generalizations of 
the ideal gas Fock-space structure to non-trivial statistics. 



5. Conclusion 

I have attempted to present, from my perspective, the main results associated with the 
still-unfolding properties of the Haldane-Shastry spin chain model, stressing the simplest 
SU{2) or S = 1/2 model. I have clearly omitted many aspects, and there is clearly much 
more yet to emerge. 

A direct algebraic treatment in the rational limit, when the "quantum group" becomes 
classical, would be particularly desirable. It can be no accident that this is the limit in 
which all the explicit results are obtained, but (disappointingly) so far a direct use of its 
algebraic properties such as the infinite 5/2+ symmetry not yet been made. An algebraic 
construction of the highest weight eigenstates (fully spin-polarized spinon gas states) 
in terms of particle creation operators ( "vertex operators" ) acting on the vacuum is also 
needed; this would closely parallel a similar treatment needed for the Calogero-Sutherland 
model. 

Other possible lines of investigation are whether a "quantum deformation" 5/2+ — ^ 
Y{sl2) can be used to calculate the hyperbolic model form factors, and what is the origin 
of integrability of the elliptic models. There is clearly much more work to be done! 
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